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Abstract 

Point sources in complex spacetime, which generate acoustic 
and electromagnetic pulsed-beam wavelets, are rigorously de- 
fined and computed with a view toward their realization. 

1. Introduction 

Acoustic and EM waves are natural candidates for multiscale 
analysis: the wave/Maxwell equations have no a priori scale. 
This motivated the idea of physical wavelets [K94a, K02a], 
localized solutions W z giving frames in solution spaces. 
z = x + iy is a complex spacetime point generalizing the 
time and scale parameters in ID wavelets. 
W z propagates and can't have bounded spatial support. 
But its source does, giving z a simple interpretation: 
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• x — (x,t) is the source center in space and time 

• y — (y, u) gives its extension around x: 

o y gives radius and orientation of launching disk 
o u gives pulse duration 

• (x, y) ~ 4D version of time & scale 

• Stability =>- y is timelike =>- z E T = causal tube. 
Thus W z is a pulsed beam (PB) with origin, direction, 
sharpness and duration controlled by y. 

Folklore: In some sense, 

W z ~ wave emitted by a complex sourcepoint z. 

Since the 1970s, complex-source beams have been used in 
engineering; see [HF01] for a comprehensive review. 
Related ideas have circulated in relativity since the 1960s, e.g., 
the Kerr-Newman solution [N65] representing charged, spin- 
ning black holes; see also [KOla, N02]. 

My motivation came from a long-time project combining QM 
with relativity via coherent-state representations, with T as 
extended phase space [K77, K78, K87, K90]. 
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• PB wavelet analysis of waves - local alternative to Fourier. 

• Associated complex sourcepoint analysis of sources. 

• Applications to radar and communications: [K96, K97, KOI]. 
Can platforms he made to launch and detect W z 's? 

Must understand the nature of the sources, making precise 
the notion of a wave emitted "from" x + iy. See also [HLKOO]. 

2. Point sources in complex space 

First define a point source in complex space C n ,n > 3. 
(n = 2 is special and must be treated separately.) 
The point source at x = G W" 1 can be defined as the source 
of the Newtonian potential G n : 



5 n (x) = A n G n (x), G n (x) = — 



where oo n is the area of the unit sphere in M n and 





is the Euclidean distance. 



Define the complex-distance function 



r{x + iy) = yj{x + iy) 2 = \[7 2 ~ 



— a 2 + lix ■ y, a = \y\ > 0. 
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Fixing y ^ 0, the branch points form an (n — 2)-sphere 

B(y) = {xeW 1 :r = a, x ■ y = 0}. 
The simplest branch cut for which f(x) = +r(x) is 

Re f > 0, which implies — a < Im f < a. 
The branch cut is the (n — l)-disk spanning B(y), 

V(y) = {x:r<a,x-y = 0}, dV = B. 
We now define the point source at x = —iy by 

S n (x + iy) = A n G n (x + iy), G n (z) = — 

uj n 2 — n 

where A n is the distributional Laplacian in x. 

5 n (x + iy) is a distribution supported in x G T>(y) [K00]. 

For even n, G n (z) is analytic whenever z 2 ^ 0, hence 

su PPa; $n( x + "iy) — B(y) f° r even n > 4. 
But for odd n, G n inherits the branch cut and 

su PPcc &n{x + iy) = T^{y) fc> r odd n > 3. 
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The distribution 5^{z) will be computed along with the time- 
dependent one for pulsed beams in M 3 ' 1 . 

3. Complex spacetime sourcepoints 
Begin with Euclidean M 4 and complexify: 

x E = (x,-u), y E = (y,t)eR 4 
z = x E + iy E = (x + iy, —u + it) e C 4 . 

Rewrite this as a complex Minkowski vector: 

z = x + iy = (z, it), z = x + iy, t = t + iu 
x = (x,it), y = (y,iu) e R 3 ' 1 

2 2,2 2 222 22 

x — x — t , y = y — u , z — z — t . 
Consider the Newtonian potential in IR 4 and its extension, 

G±(x E ) = — - — 5—5-, G^(z) = — - — t— -. 

Gi(x E ) is a fundamental solution for the Laplacian, 

/\±G 4 {x, u) = (Aa, + dl)G 4 (x, u) = 5 A {x, u). (£) 
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Naive question: Does u — > it gives a fundamental solution for 
the wave operator, i.e., 

DG A (x, it) = (A x - d 2 ) G 4 (x, it) ~ 5(x, t) ?? 
Answer: No! 

• G 4 is singular on x 2 = 0, so it must be defined in M 3,1 . 

• When properly defined, it will be sourceless. 

• Causality makes sense in R 3 ' 1 but not in M 4 . 
To find the correct definition, note that 

-z 2 = t 2 - z 2 = t 2 - f 2 = (r + f)(r - f), 

which gives the partial-fractions decomposition 



G 4 (z)=G + (z)-G-(z), G±{z) 1 



8n 2 r t =F r 



Since formal differentiation gives 

OG ± (z) = 0, 
G ± (z) can have sources only at x e V or r = ±f. 



Define the time-dependent radiation pattern 

2ttR^ 



r=Fr (t TP) + i(u =F q) 



r = p + iq, p > 0, —a < q < a 



Now (p, q) are oblate spheroidal coordinates in R 3 whose 
level surfaces are £>-confocal ellipsoids and hyperboloids. 



In cylindrical coordinates with X3 = y ■ x and p = a/v 2 — x 
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E p = {constant p > 0} = < x : 



x% 

+ -1 = 1 



a 2 + p 2 p 2 



H q = {constant q 2 < a 2 } = |a; : 



p 2 4 _ x 

a 2 — q 2 q 2 



EpS are wave fronts, and HqS give orthogonal flow. 
In the far zone r ^> a, they become spheres and cones: 

f = \/r 2 — 0? + liar cos 6 ~ r + ia cos 
-^p i r = #<z {cos6> = q/a] 



2ttR 



t=±r 



(t =f r) + i(u =f acos6*) -u=Facos#' 



(*) 



(I) y is timelike: |u| > a, is a smooth pulse around ±y 
peaking at t — ±r, with duration T(x) = \u =F acos#| and 
elliptical radiation pattern of eccentricity a/\u\. 



(II) y is spacelike: \u\ < a, G ± is singular on the cone 
cos 9 = ±u/a at t = ±r, with a hyperbolic radiation pattern. 

(III) y is lightlike: \u\ = a, G ± is singular on ray cos 9 — ±1 
at i = ±r, with a parabolic radiation pattern. 

Only (I) gives a reasonable PB, though (II, III) should be of 
interest otherwise since G\{z) is holomorphic for all z 2 ^ 0. 
We therefore assume 

y G V± = {(y,iu) : ±u > \y\} = future/past cone, 

so that z belongs to the causal tube T = T + U 71, where 



71 are famous in physics (quantum field theory, twistors) as 
well as mathematics (Lie groups, harmonic analysis). 
Relation to propagators: Fix y e V ± and write 

G K (x±iO)= lim G K (a;±^), « = +, -. 

Then (4) gives (with P= principal value) 



71 = {x + iy : y e V±} = future/past tube. 



G ± (x - iO) 



G ± 



(x + iO) 



8n 2 r t =F 



87r 2 r t =F r 




+ 



2<5(t =F r) 

87rr 
i<5(t =F r) 

8nr 



8 



Huygens' principle requires the Minkowskian limit 

G£(x) = G ± {x- iO) -G ± (x + iO) = %5 ^ r \ ( M ) 

These are the retarded and advanced propagators, fun- 
damental solutions vanishng for ±t < 0: 

□G£(s) = -i5(x)5(t) = 5(x), (TL) 

where the last equality comes from our volume element, 

x = (x,it) e K 3 ' 1 =>- dx = idxdt. 

• (H) is the desired hyperbolic counterpart of (£). 

• G M = G+j — G M is the Riemann function: 

Gu (x) = ^ (t ~ r) _ t6 ( t + r ) 
Anr Airr 

□G M (x) = 0, G M (x,0)=0, d it G M (x,0)=5(x). 

• Causality comes with a choice of branch and has no 

meaning for G 4 (z): f — > — f =>- G ± — > — G T . 

(M) is typical of hyperfunction theory, representing distri- 
butions as limits of differences of local holomorphic functions 
(in general, sheaf cohomology classes) [K88]. 
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We now define the point source at x = —iy as 



5(z) = DG ± (z), z = x + iy. 



(*) 



• This is identical for G + and G ~ : 

z 2 = =>- x 2 = y 2 < and x ■ y — 0, 

but x • y 7^ since both vectors are timelike 
Hence G4 is holomorphic in T and 

UG + (z) - UG-(z) = DG 4 {z) =0. 

• By (Ti), the Minkowskian limit of 5 is 5: 

<5 M (x) = 5(x - iO) - 5(x + iO) = 5(x). (M 1 ) 

• Since r =)= f 7^ in T, G ± is singular only when x EV and 

supp 5(x + iy) = {(a;, it) : x £ V{y)} = V(y) "iy e V ± . 



V{y) is the world tube swept out in 1 ' by T>{y) at rest. 
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4. Computing the source distributions 

Consider the function on T defined by 

W{z,iT) = 9< " T ~f\ K = ± (W) 

47rr 

with g(Y) an analytic signal, holomorphic for Im r ^ 0. 
We want to compute the source distribution 

s(z) = n x w(z). 

Note that 

g(r) = JL =► = =► = ~5{z) 

g {r) = -1 W^) = G 3 (z) = 5(z). 

• Im (r ± f) 7^ in T =>- supp 5 C V{y). 

• W is doubly singular on X>, where g(r — nf) has a jump. 
Regularize W with the Heaviside function 6. Let e > and 

W(z) outside E £ 
inside E £ . 
Then the regularized source, defined by 

s e (z) = nw e (z), 



W £ (z) = e(p-e)W(z) 
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is supported on the world tube E £ C M 3 ' 1 of E £ . 
A computation gives 

Jp 2 + a 2 )g Je 2 + a 2 )(rg p - g) eg 



\f\ z S £ = 5' 



Aur 



+ 5- 



2nP 



+ 5 



2iir 



where 



S(p - e) = 9'(p - e), g p = d p g = -ng'{r - «f). 
S £ acts on a test function f(x) (no smearing needed in t) by 
(S £ ,f) = J dxS £ (x + iy,ir)f(x) 



p 2 + q 2 

with dx = dp dq d<p. (Vol) 

a 



Integrating by parts in p and simplifying gives 



(Se,f) = 



e 2 + a 2 ra 



2a 



dq 



9 P f gf gf P 



p=e 



where f(p, q) is the mean of f(p, q, 0) over 0. 

But g p = —ig q as g is holomorphic. Integrating by parts in q, 



'gf 


r=a 

aa f a 7 
dq 


gff 






r 


. _ a J-a 

r=a 


f 



(0) 
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with a = e + ia and a the north and south poles of E £ , 

/(f) = f(p, q) (no analyticity implied in p + iq) and 

ff(r) = dff(f) = \i.d p - id q )f(p, q). (f f ) 

Taking the limit e — > now gives the action of S: 

(S,f) = -g(r, a)f(0) + 2iaf^ g(r, q)f,(iq) (9) 

Jo 1 P= e 

where 

<K r > ?) = \ [g( r - n) + g( T + 

is the average of #(t — nr) over the jump at f = ±iq e D and 
we used the continuity of / and its derivatives across V, 

/(-»?) = /(*?), /?(-»?) = -/f(*9) 
/(±ia) = /(0), /r(0) = 0, 

which ensures that the integral in ("v 1 ) is defined. 

As a check, note that ("v 1 ) gives the correct value as y — > 0: 

(5,/>--</(r)/(0) =► 5(a;,<r) = -<7(T)(y(aj). 
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5. Interpretation 

Let us examine the meaning of the right side in 



_ f a d gff 

a ./_„ f 



(0) 



• The boundary term is a pair of real sourcepoints at the 

poles f — a, a (p — 0, £3 = ±e) of -E £ modulated by 

• Since d p is outward-orthogonal to E £ , the term with <9 p / is a 
double layer on E £ , modulated by g; see Eq. (f f ). 

• 9 9 is tangent to E e , representing a flow from a to a. Hence 
the term with d q f is a flow on modulated by g. 

• Since g{r) is an analytic signal [K94], the modulation by 

g = g (j T f) = g(t =f e + i(u =F g)) 

is strongly amplified if ±qu > and diminished if ±gw < 0. 
This makes W(z) a unidirectional beam. 

Taking into account (Vol) gives the unsmeared form 



S e (z) = W{z) 



iS(f — a) — iS(f — a) — 2 



f 



5(p — e) 9r 



where 



5(r — e ± ia) = 5(p — e)5(q ± a). 
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Similarly, ( < s?) gives the unsmeared form 
S(z) = -g(T,a)5(x) -g(r,q) 



a 2 5(p) 
2niq 3 



df 



In cylindrical coordinates, this becomes 



S(z) = -~g(T,0)6(x)-~g(T,p) 



Q(q - p)5(x 3 ) /a 
2n V V p 1 V' 



d p - id 3 



where 



1 



a 2 — p 2 ) 



Letting g = — 1 gives the static complex sourcepoint 



5(z) = 5(x) + 



0(a — p)5(rr 3 ) / a 



2ir y /a 2 -p 2 VP 



which is simpler than the form derived in [K00], being local 
while the latter had a subtraction. 

Letting g{r) = 1/2ht gives the complex spacetime point source, 



S(z) = 



tS(x) r6(a — p)S(x 3 ) fa 

+ 2lXZ 2 ^^ V " 



where 



z = (z,ir) e V 



2 2 2 / n 
P —a — T f: 0. 
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• The above method works in R™' for all odd n > 3. 

• An interesting connection has been obtained between the 
Euclidean source distribution 5 n+ i and solutions of the homo- 
geneous wave equation in R™' 1 , based on the fact that the 
latter are spherical means over (n— l)-spheres and <5 n+1 is 
supported precisely on these spheres for odd n [K00]. 

• These results extend to even n > 2 by a variant of Hadamard's 
method of descent relating 5 n to 5 n+ i and to 5 n +i t \. 

The cases 62 and 8^1 are special since G2 is logarithmic, but 
the results are similar. 

6. Making waves 

Acoustic wavelets are generated from a 'mother' PB W ± by 
W±{x) = W ± (x' - z) = W ± {x -x + iy). 

For suitable g, subfamilies of W^s are frames, giving a pulsed- 
beam analysis-synthesis scheme for general waves as a local 
alternative to Fourier analysis. 

In [K94], the W(z)'s are (fractional) time derivatives of G ± (z). 
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Electromagnetic wavelets are frames of dyadic PB propa- 
gators Gz(x') = G ± (x'-z) constructed as follows [K94, K02]. 

1. An electric & magnetic (E&M) dipole with dipole moments 
p e ,p m at x — —iy gives an E&M polarization density 

P(z) = P e (z) + iP m (z) = p S(z) , p = p e + ip m e C 3 . 

2. This gives retarded/advanced E&M PB Hertz potentials 

Z ± {z) = Z±{z)+iZt n {z)=pG ± {z) 
UZ ± (z) =pUG ± (z) = P ± {z). 

3. These in turn generate EM fields F ± = E ± + iB ± given by 

F ± (z) = V x V x Z ± {z) - d it V x Z ± (z). 

4. The dyadic 'mother' PBs are now defined by 

G ± (z)p = F ± {z). 

• All the Gelds are real! The self-dual pairings E + iB, etc. 
merely exhibit the (complex!) structure of EM. 

• Gives local analysis/synthesis schemes for fields and sources. 

• PB analysis seems more natural for EM than scalar waves 
since the real and imaginary parts have direct interpretations. 
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Applications to radar and communications have been proposed 
[K96, K97, KOI] whose utility would be greatly enhanced if 
sources can be constructed to launch and detect EM wavelets. 
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